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Abstract 

We investigate the critical temperature of an interacting Bose gas confined 
in a trap described by a generic isotropic power-law potential. We compare 
the results with respect to the non-interacting case. In particular, we derive 
an analytical formula for the shift of the critical temperature holding to first 



order in the scattering length. We show that this shift scales as iV 3 ( n + 2 ) , where 
iV is the number of Bosons and n is the exponent of the power-law potential. 
Moreover, the sign of the shift critically depends on the power-law exponent 
n. Finally, we find that the shift of the critical temperature due to finite-size 

_ 2n 

effects vanishes as ./V 3 (™+ 2 ) . 
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Nowadays more than twenty experimental groups have achieved Bose-Einstein condensa- 
tion (BEC) in trapped alkali-metal atoms by using different geometries of the confining trap 
and atomic species. 1-3 The experimental results on these dilute gases of weakly- interacting 
Bosons are in quite good agreement with the theoretical predictions both above and below 
the BEC critical temperature (for a recent review see Ref. 4). 

An interesting problem is the behavior of the critical temperature T c as a function of 
the number N of Bosons. Unfortunately, the experimental data of the critical tempera- 
ture are not very precise but they seem to qualitatively confirm the predictions. In many 
experiments with alkali-metal atoms the trap can be described by a harmonic potential 
U(r) = (l/2)mu 2 r 2 . In this case, it is well known that in the large N limit the critical tem- 
perature T^> of an ideal gas scales as iVi, 5 while the shifts of the critical temperature due to 
the finite number of Bosons and the interaction scale as and iVs, respectively 6,7 Note 
that the shape of the trapping potential plays a decisive role for the critical temperature, as 
we have recently shown with a dilute gas of hydrogen atoms in a Ioffe trap. 8 

An important class of trapping potentials is provided by power-law potentials U(r) = 
A r n . Power-law potentials have been proposed to cool the Bose gas in a reversible way 
by adiabatically changing the shape of the trap at a rate slow compared to the internal 
equilibration rate. 9-11 For these potentials it has been shown that the critical temperature 
T^> of an ideal Bose gas scales as iV 3 (™+ 2 ) in a 3-dimensional space 5 and as N d ( n + 2 ) in a 
generic d-dimensional space. 12 

In this paper we extend previous results of an ideal Bose gas in power-law potentials. 
In particular, we study the role of the inter-atomic interaction and finite-size effects. We 
derive analytical formulas for the total number of particles and the critical temperature of 
BEC phase transition. 

The starting point of our investigation is the semiclassical density profile p^°^(r) of an 
ideal Bosons in a trapping external potential U(r). In the grand canonical ensemble, the 
density profile is given by 
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P (0) (r)=(^) S Ss/2(e^»/ fcr ) , (1) 

where T the absolute temperature with k the Boltzmann constant, \x is the chemical potential 
and 

oo i 

B r {z) = Y,-, (2) 
i=i 1 

is the so-called Bose function. This result is the generalization of the formula for an 
ideal Bose gas in a box of volume V . It shows that, in the semiclassical limit, the non- 
homogeneous formula is obtained with the substitution \i — > fj, — £/(r), also called local 
density approximation. 7,12 

In the case of weakly-interacting Bosons, the previous formula must be substituted by 
the following one 

p(r) = (^) i S3/ 2 (e^ + ^ fcr ) , (3) 

where g = Anti^a/m is the interaction coupling constant fixed by the s-wave scattering 
length a. This equation is a mean-field self-consistent equation that can be derived, within 
the semiclassical approximation and above the critical temperature T c , from the Bogolibov- 
Popov equations, which describe the elementary excitations of the Bose gas. 4 ' 7 The chemical 
potential /x is fixed by imposing the following normalization condition 

N = J p(r) d 3 r , (4) 

where N is the total number of Bosons. Note that below the critical temperature T c , the 
density of the left-hand side in Eq. (3) describes only the non-condensed thermal cloud, while 
the density in the right-hand side is the total density, namely condensed and non-condensed 
density. Nevertheless, to calculate the BEC transition temperature T c it is sufficient the 
study the non-condensed thermal cloud. In fact, by imposing for the density profile (3) the 
normalization condition (4) one finds the total number of particles N as a function of T and 
/i. At the critical temperature T c , the chemical potential fi is minimum (/i = in the large 
N limit) and thus one gets the critical temperature T c as an implicit function of N. 4 ~ 7,12 
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Let us consider a generic isotropic power-law potential given by 



(5) 



where n is the power-law exponent and A is the trap constant. Here we have introduced an 
energy parameter hu and a length parameter r , which can be chosen as r = {%/ {muoo)) 1 ^ 2 . 
By Taylor expanding the Eq. (3) to the first order in the coupling constant g one finds 



p(r)=p(°)(r)- 



2g ( mkT s 



B 1/2 (< 



(VL-U(t))/KT 



)p (0) W, 



(6) 



kT \2irh 2 / 

where p^(r) is the non-interacting density given by Eq. (1). Then, by imposing the 



normalization condition one gets 

3(n + 2) o 

V^o/ 2tr(f + 1) ^ v 

where T(r) is the Euler factorial function and 



l-2s 



(*T)* 



3(n + 2) (e^ T ) 
2n 



, (7) 



.3 .i 



(8) 



lJ = 1 i2j*(i + j) r 

In Eq. (7) the total number iV of Bosons is a function of temperature T, chemical potential 
/i and interaction constant g, namely N = N(T,fx,g). Setting g = we have the number 
of particles for an ideal Bose gas. In this case, to find the critical temperature we put 
fx — and obtain 



&T C (0) = c(n) huj iVrai , 



(9) 



where c(n) is a numerical coefficient, 



c(n) 



2§r(| + i) 



2n 
3(n+2) 



(10) 



.2Sr(J + i)c(J + i; 

and £(r) = B r (l) is the Riemann ^-function. As previously stated, this formula was obtained 
for the first time by Bagnato, Pritchard and Kleppner. 5 

To determine the effect of the inter-atomic interaction on the critical temperature, we 
expand the Eq. (7) around fx = 0, T = T c (0) and g = 0: 

N(T, fx, g) = N(Tf\ 0, 0) + 0, 0) (T - 7?>)+ 



+ ^(T?\0,0)l*+^(T?\0,0)g. (11) 

At the critical temperature (T = T c ) of the interacting system, for large N, the chemical 
potential can be written fx = 2gp^(0) as suggested by Eq. (3). 4,7 In this way, from Eq. (1), 
Eq. (9) and Eq. (11) we get the main result of the paper 

= ~d(n) (£) AT^T , (12) 

where 

2lr(| + i) 



M , 2l(c(f + |)C(l)-^(f 
{n) = **C(J + §)(* + §) 



2-r(J + i)c(| + l) 



3(n+2) 



(13) 



and S'(r) = «SV(1). An inspection of Eq. (12) shows that the shift due to the interaction can 
be either negative or positive, depending on the sign of the scattering length a and of the 
coefficient d(n). As shown in Tab. 1, d(n) diverges for n = 6 and goes to zero as n — > oo. 

n 

The temperature shift (12) scales as jV 3 (™+ 2 ), where n is the exponent of the power-law 
potential. It means that increasing n one can amplify the shift of the critical temperature 
due to the interaction. This effect is shown in Fig. 1, where we plot the critical temperature 
T c as a function of the number N of bosons for n = 2 and n — 4. It is important to observe 
that in our approach the shift is always associated with a change in the central density 
produced by inter-atomic forces. The relationship between the critical temperature and the 
corresponding value of the density in the center of the trap is unaffected by the interaction. 
Moreover, the shift vanishes for a homogeneous system where the density is uniform. Note 
that with n = 2 one recovers the shift-formula obtained by Giorgini, Pitaevskii and Stringari 
for a dilute Bose gas in harmonic potential. 7 

Another interesting aspect of BEC in trapping potentials are finite-size effects. In an ideal 
gas, BEC starts at the temperature for which the chemical potential fx reaches the energy 
of the lowest solution of the Schrodinger equation. In the case of a power-law potential 
this energy is a n hu , where a n is a coefficient depending on the power-law exponent n. For 
example, in the case of the harmonic potential we have a.2 = 3/2. In general, the coefficient 



a n must be numerically evaluated and only in the large N limit one can set [i — at the 
transition temperature. We can estimate finite-size effects on the basis of Eq. (11) by setting 
g — 0, ii — a n htuo and T = T^ FS \ In this way we obtain 



(o) 



= -e(n)a n N 



(14) 



where 



e ( n ) = 773 i 3\/ 3 r 3 



C(f + I)(f + 1) 



2§r(| + i) 



2n 
"3(n + 2) 



(15) 



.2*r(J + i)c(i + i; 

Numerical values of the coefficient e(n) are reported in Tab. 1. The temperature shift 
(14) originating from finite size effects changes sign as at n = 6, where the coefficient e(n) 
diverges, and, as expected, the shift vanishes in the large N limit. In particular, with a 
fixed number iV of Bosons, one can strongly reduce this finite-size shift by increasing the 
power-law exponent n as shown in Fig. 1, where we have set a 4 = a 2 = 3/2. The finite-size 
shift predicted by Grossmann and Holthaus 6 for a Bose gas in harmonic potential is obtained 
from Eq. (14) with n = 2. 

In conclusion, we have derived two remarkable formulas for the critical temperature 
of a weakly-interacting Bose gas in a generic isotropic power-law potential. One formula 
gives, to the first order in the scattering length, the shift of the critical temperature due 
to the inter-atomic interaction. The temperature shift is strongly enhanced by a large the 
power-law exponent but the sign of the shift depends in a non trivial way on the power-law 
exponent. The other formula shows that, as expected, the finite-size effects, namely setting 
the chemical potential equal to zero at the Bose-Einstein transition, are negligible in the 
large iV limit. Moreover, with a fixed number N of bosons, the finite-size effects are reduced 
by a large power-law exponent. 
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1 2 3 4 5 6 7 8 


c(n) 
d{n) 
e(n) 

2n 
3(n+2) 

n 

3(n+2) 


0.56 0.94 1.14 1.25 1.31 1.35 1.04 1.00 
0.64 1.33 -0.23 2.41 7.24 oo -41.57 -42.96 
0.40 0.50 0.29 1.11 2.45 oo -2.96 -1.61 
2/9 1/3 2/5 4/9 10/21 1/2 14/27 8/15 
1/9 1/6 1/5 2/9 5/21 1/4 7/27 4/15 



TAB. 1. Numerical coefficients as a function of the power-law exponent n. In the limit 

n -»• oo one gets: c(n) -»• 1.27, d{n) -> 0, e(n) -> -0.29, 2n/(3(n + 2)) -»• 2/3, and 
n/(3(n + 2)) -»• 1/3. 
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FIG. 1. Critical temperature T c as a function of the number N of Bosons. Power-law exponent 
n = 2 on left and n = 4 on the right. Ideal Bose gas (dashed line), interacting Bose gas (long dashed 
line), interacting Bose gas with finite-size corrections (full line). Interaction strength a/ro = 10 -2 . 
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